High mobility topological interface state probed by terahertz measurements 
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We report gated terahertz cyclotron resonance measurements on epitaxial Bi2Se3 thin films capped 
with In2Se3. The properties of a single topological interface state are measured from the vicinity of 
the Dirac point to above the conduction band edge. A precipitous drop in the scattering rate with 
Fermi energy indicates decoupling of the surface states from bulk states and a shift of the Dirac 
cone towards mid-gap. Near the Dirac point, the mobility is 3500 cm 2 /V-s with measured potential 
fluctuations of 60 meV. 
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The topological insulator is a unique state of mat- 
ter that is a bulk insulator possessing a topologically 
protected metallic surface state (TSS) of massless par- 
ticles known as Dirac fcrmions. Key properties of this 
spin textured TSS are suppression of back scattering and 
an intrinsic magneto-electric effect* 1 ^ The experimental 
realization of this state of matter in Bi 2 Se 3 and other 
compounds has sparked great interest owing not only 
to their potential use in spintronic devices but in the 
investigation of the fundamental nature of topologically 
nontrivial quantum matter Bi2Se3 is one of the most 
important topological insulators due to a single Dirac 
point existing above the valence band edge and a large 
bulk bandgap of 300 meV offering the prospect of room 
temperature applications^ Unfortunately, unintentional 
doping of these materials has led to large bulk densi- 
ties at the surface that mask the electronic properties of 
the surface state.— ~— Capping materials provide a way 
to control this doping, both chemically and as a plat- 
form for gating, but the usual surface probes that have 
studied the vacuum/TSS interface stat o 1 ' 2 ' 12 ^ 3 can not 
measure buried interfaces where critical Dirac cone prop- 
erties are expected to change d 4 ' 15 In this letter we report 
gated terahertz (THz) cyclotron resonance measurements 
on passivated Bi2Se3 thin films that directly measure the 
scattering rate and cyclotron mass of a single topological 
surface state from above the conduction band edge to the 
vicinity of the Dirac point. 

Transmission measurements through a thin film of 
topological insulator probe all carriers, which can include 
topological surface states, accumulation layers, and bulk 
carriers as depicted in Figure QJa). Modulating the topo- 
logically protected interface state of a capped him with a 
scmitransparent top gate while measuring the modulated 
terahertz transmission allows extraction of its properties, 
even in the presence of much larger background conduc- 
tivity contributions. Such measurements performed in 
magnetic held give not only the carrier density and scat- 



tering rate but also the sign of the charge and cyclotron 
mass of the carriers. For negative gate biases that move 
the Fermi level of the surface state well below the con- 
duction band edge, Dirac cone electrons are expected to 
exhibit a small cyclotron mass m c = HRf/vf since the 
Fermi wavevector kp goes to zero while the Fermi veloc- 
ity vf remains finite at the Dirac point 

An epitaxial Bi2Se3 60 quintuple layer film and a 10 
nm In2Se3 capping layer were deposited without break- 
ing vacuum onto a sapphire substrate* 1 ^ The dielectric 
parylene-C and a NiCr semitransparent gate were sub- 
sequently deposited (see the Supplemental Materials for 
details 2 ^). The device is depicted in Figure D^b). 

Zero-gate normally incident transmission measure- 
ments are shown in Figure [IJc,d). Transmission T 
through a thin film is related to the two dimensional 
conductivity a by T ~ \n su b + 1 + Zq<j\~ 2 , where n su & 
is the index of refraction of the substrate and Zq is 
the impedance of free space, a is a sum of all con- 
tributing channels. The Drude conductivity in a mag- 
netic field expressed in the circular polarization basis is 
a = (u)p S / Zq) / — i(us ± cl> c )), where the two dimen- 
sional plasma frequency is u ps = Z^ne /m and u> is the 
radiation frequency. The cyclotron frequency ui c is re- 
lated to the carrier charge e and cyclotron mass m c by 
lu c = eB I (cto c ), where B is the applied magnetic field. In 
zero magnetic field, this conductivity expression depends 
only on the two dimensional carrier density n, effective 
mass to, and scattering rate 7. 

The Fourier transform transmission measurement of 
Figurc[]Jc) was performed in zero magnetic field with un- 
polarizcd incident light. The cyclotron resonance trans- 
mission measurement is shown in the inset of Figure [He). 
It was performed with circularly polarized incident light 
at a fixed laser frequency as a function of magnetic field 
applied normal to the film. The cyclotron resonance mea- 
surement exhibits a single Lorentzian absorption. The 
resonant position located in positive magnetic field cor- 
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responds to n-type carriers of mass m c « 0.19 mo- The 
Faraday angle Op reported in Figure [TJd) was measured 
with polarization-modulated incident light at a fixed fre- 
quency, a method detailed elsewhere^ Fits to the zero- 
gated transmission data in Figure [ljc,d) using a single- 
fluid Drude model for the conductivity, depicted by the 
blue curves, reveal a large carrier density n « 10 13 cm~ 2 
and a cyclotron mass m c > 0.19 mo that is much larger 
than the conduction band edge mass m c = 0.15 mo^^ 
Deviations between the Op data and fit are outside of 
measurement error, hinting at small contributions from 
other conducting channels. 

Transport measurements in similarly grown epitax- 
ial films, performed over a wide range of thicknesses, 
show the bulk contribution to the conductivity is a small 
fraction compared to the surface layer contributions.— 
Therefore, the observed large carrier density is likely as- 
sociated with a top and/or bottom surface accumulation 
layer due to defects. 

Modulating the top gate voltage at ~ 1 Hz and mea- 
suring the difference in cyclotron resonance (A-CR) re- 
veals the top topological interface state properties. The 
largest negative average gate voltage represented by the 
red curve of Figure dja) shows a resonant component 
at B w 2 T. This corresponds to half the conduction 
band edge mass i 6 ' 7 ' 22 Such a small mass is evidence of 
the topological interface state Dirac cone and a Fermi 
level of approximately 50 mcV above the Dirac point, far 
below the conduction band edge of 190 meV observed on 
freshly cleaved Bi 2 Se3 in vacuum . 16 ! 17 

The small mass associated with the top surface estab- 
lishes that the large accumulation layer dominating the 
zero-gate response is on the bottom surface. Changing 
the gate from zero to the maximum of -170 V transfers 
4.4 x 10 12 e/cm 2 , a fraction of the 10 13 e/cm 2 associated 
with the zero-gate measurement. Gating the top topo- 
logical surface state to the vicinity of the Dirac point 
requires removing any accumulated charge in addition 
to depleting the surface state and a region of the bulk 
near the surface. Given the gate limitation, this is not 
possible if the top interface is accumulated. The large 
accumulation layer at the bottom surface includes a bot- 
tom surface topological surface state. These two large 
mass contributions are indistinguishable in our differen- 
tial measurements since the bottom surface is not signif- 
icantly modulated by a top gate. 

The observed mass dependence on gate voltage is con- 
sistent with a Dirac cone surface state below the con- 
duction band edgei 16 ! 17 In order to gain further insight 
into the gating behavior of the surface state, the optical 
response of the film is modeled in terms of the surface 
state dispersion as measured by angular-resolved pho- 
toemission spectroscopy (ARPES) i 16 ' 17 The model in- 
corporates three n-type conductivity contributions dis- 
tributed in the film (details are provided in the Supple- 
mental Materials 2 ^) : one term for the top topological sur- 
face state (TSS), another for the top accumulation layer 
and/or interim bulk carriers (IB), and the third term 




FIG. 1: Zero-gate normal incident transmission mea- 
surements on epitaxially grown B2Se3 films (a) A band 
bending schematic showing the evolution of the conduction 
band edge relative to the bulk chemical potential (green), and 
the corresponding Dirac cone Fermi level associated with the 
vacuum interface, when a variable gate voltage is applied to 
the top surface. Purple through red correspond to increas- 
ingly negative gate voltages. For large negative voltages, a 
depleted region exists decoupling the surface state from the 
bulk (red). The bottom surface (black) is accumulated due to 
defects, (b) A schematic represents the measured device, (c) 
Measured 6 K Fourier transform (FTIR) transmission spec- 
trum (black). A single- fluid Drude model fit (blue) gives 
7 = 7.3 THz and lj ps /(2it) = 27.9 THz that translates to a 
carrier density n = 1.23 x 10 13 cm -2 for a mass of m = 0.21. 
Oscillations in the FTIR data are due to Fabry-Perot etalon in 
the sapphire substrate, (inset) Measured 10 K cyclotron reso- 
nance transmission at o;/(27r)=0.74 THz. (d) Complex Fara- 
day angle measured at 10 K and u>/(27r)=0.74 THz (black). 
A single-fluid Drude model fit (blue) gives 7 = 7.4 THz, 
n = 1.1 X 10 13 cm -2 , and m c = 0.21 . The red curves in 
figures (c) and (d) are fits using a three-fluid Drude model 
described in the text. 



for the bottom surface accumulation layer (BSAL). Each 
conductivity contribution is expressed within a Drude 
model with parameters for carrier density n, scattering 
rate 7, and cyclotron mass m c . The BSAL is not mod- 
ulated by the gate and its Drude parameters are cho- 
sen to be consistent with the zero-gate measurements. 
To reduce the number of parameters associated with the 
gate modulated carriers, the amount of charge a gate 
transfers on or off the topological surface state is deter- 
mined by solving for the band bending in the film us- 
ing Thomas-Fermi screening theory. The IB parameters 
are gate-dependent where m c and 7 are free parameters, 
but the change in carrier density is determined by the 
band bending model. The gate dependent TSS parame- 
ters n and m c are set by the ARPES measured dispersion 
relatio n 16 ' 17 while 7 is the only free Drude parameter. 
The differential cyclotron resonance measurements are 
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FIG. 2: Gated magneto-optical transmission measurements at 10K: Average gate values V g for each plot are coded 
according to the colored bar at the top of the figure. All optical signals are differences between two gate values, V g ± AV/2, 
measured at every integer B-field value up to ±8 T. Shown on each data graph is the radiation frequency u)/(2n) and AV 
value, (a) Differential cyclotron resonance data normalized to the average with a standard deviation of the mean of 5 x 1CP 4 . 
(b,c) Differential Faraday angle data measured with a standard deviation of the mean of 30 ^rad (d-f) The modeled optical 
responses are detailed in text. The colors of the modeled curves correspond to the same V g and AV' values appearing in the 
respective data sets (a-c). 



mainly sensitive to the absorptive (real) part of the con- 
ductivity. Additional information is gained from the dif- 
ferential complex Faraday angle data reported in Figures 
[U(b,c). The real (imaginary) part of 9p is related to the 
reactive (absorptive) part of the conductivity, 2 * 

Figures UK d-f) show the modeled differential optical 
signals. The red curves in Figures [ljc,d) show the zero- 
gate Fourier transform spectroscopy and Faraday angle 
model results. A comparisons to other A-CR measure- 
ments at higher frequencies as well as gated spectroscopy 
data are shown in the Supplemental Materials^ Figure[3] 
summarizes the Drude parameters which reproduce the 
data sets. All the features are reproduced over a wide 
range of applied voltages, magnetic fields, and frequen- 
cies. 

In the A-CR data of Figure [2]( a), the observed peak- 
dip-hump structure evolution at negative gate voltage is 
reproduced by the model, showing that the cyclotron 
mass m c = hkp/vp is consistent with the Dirac dis- 
persion measured by ARPES (a more thorough intuitive 
discussion of the model is presented in the Supplemental 
Materials 2 ^). However, the modeled differential response 
at the highest negative gate voltages is larger across all 
data sets. The differential optical signals are proportional 
to changes in the optical spectral weight, ne 2 /m c . The 
Dirac cone dispersion predicts n/m c oc kpvp which goes 
to zero at the Dirac point as do both n and m c . In the 



presence of potential fluctuations^^— these quantities 
saturate to some non-zero rms values as the Dirac point 
is approached. This saturation of the spectral weight in 
the vicinity of the Dirac point reduces the differential 
spectral weight below that predicted by a Dirac cone dis- 
persion. The observed smaller differential optical data 
near the Dirac point demonstrate deviations from Dirac 
cone expectations, consistent with potential fluctuations. 

These deviations provide a means to characterize the 
potential fluctuation energy. The first significant devia- 
tion between model and data begins with average voltage 
V g = —105 V with its associated lowest voltage of —130 
V. From Figure E](b), the TSS Fermi level corresponds to 
~ 60 meV above the Dirac point. This estimate of the 
potential fluctuations agrees well with other estimates, 25 

The TSS scattering rate, shown in Figure E](e), mono- 
tonically decreases with Fermi level transitioning to a 
much flatter behavior at ~ 70 meV above the Dirac point, 
providing another estimate of the puddling energy scale. 
The independently measured mass and scattering rate of 
the topological surface state translate into a mobility, re- 
ported in Figures |3Jf). It is largest in the vicinity of the 
Dirac point, 3500 cm 2 /V-s. 

The scattering rate suffers a precipitous drop at ~ 70 
meV below the inferred conduction band edge that is 
based upon the ARPES dispersion. This behavior is very 
different from that observed in graphenei 24 i 27 The scat- 



4 




FIG. 3: Model parameters: Three n-type Drude terms 
model the reported differential optical measurements. The 
BSAL parameters are n — 1.0 x 10 13 cm~ 2 , m c — 0.25 mo, 
and 7 = 9.4 THz. The red through purple dots on all graphs 
depict the same average voltages applied as in Figure [2J a). 
The solid orange lines depict the depletion onset and the red 
dashed lines depict the conduction band edge, inferred from 
the ARPES measured dispersion, (a) The fraction of charge 
moved by the gate in the IB (blue) carrier contribution and 
the TSS (black), where Ept = 15 meV (screening length 
\ 3 — 13.8 nm), shown as a function of V g . (b) The TSS Fermi 
level with respect to the Dirac point inferred from ARPES 
measurements, (c) The cyclotron mass in units of the bare 
electron mass mo and (d) the scattering rate are shown for the 
IB (blue) and TSS (black) carrier contributions, (e) Scatter- 
ing rates shown as a function of the TSS Fermi level, (f) The 
mobility, calculated from the measured mass and scattering 
rate, shown versus the TSS number density. 

tering rate of the topological surface state is expected to 
decrease as it fully decouples from the bulk carriers below 
the conduction band edge, the point at which a depletion 
layer forms between the bulk and surface carriers. 

When a band- inverted semiconductor like Bi2Sc3 forms 
an interface with a trivial insulator, the Dirac point is ex- 
pected to shift relative to the conduction band edge com- 
pared with the vacuum interface studied by ARPES i 14 i 15 
The shift depends upon the size of the band gap as well 
as the relative alignment of the bands that in part de- 
pend upon the electron affinities and any potential steps 



caused by interface states^ These parameters are not 
well known for Iii2Se3 or Bi2Se3. Given the uncertainty 
in the parameters, an upward shift of ~ 70 meV in the po- 
sition of the Dirac point for the Bi2Se3/In2Se3 interface is 
reasonable (see Supplemental Materials^). Within this 
picture, the location of the step in scattering rate is natu- 
rally explained, and the width of the step is a consequence 
of potential fluctuations, ~ 30 meV. A smaller value of 
the potential fluctuations is expected near the conduc- 
tion band edge than observed in the vicinity of the Dirac 
point since there is enhanced screening by the higher den- 
sity of surface carriers as well as screening contributions 
from bulk carriers j 24 i 29 

This ability to shift the Dirac point by capping pro- 
vides a new flexibility in tailoring the topological surface 
states. In the case of Bi2Se3/In2Se3, the gap between the 
Dirac point and the valence band is increased providing 
a wider window of electron-hole symmetry. 

The gate modulated magneto-optical techniques pre- 
sented here provide a powerful means to characterize 
the topological surface state even when beneath other 
materials. The technique was used to directly measure 
the scattering rate and cyclotron mass of the buried 
Bi2Sc3/In2Se3 interface state from above the valence 
band to the vicinity of the Dirac point. A precipitous 
drop observed in the scattering rate indicates decoupling 
of the surface from bulk carriers and a shift of the Dirac 
cone by ~ 70 meV towards mid-gap. Conductance fluc- 
tuations are estimated to be ~ 30 meV in the vicinity of 
the conduction band edge and ~ 60 meV near the Dirac 
point where the mobility attains its highest value, 3500 
cm 2 /V-s. 
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SUPPLEMENTAL MATERIAL 
I. METHODS 

A Bi2Sc3 60 quintuple layer film was grown epitaxi- 
ally onto a 0.5 mm thick sapphire substrate of area lxl 
cm squared A 10 nm thick In 2 Se3 capping layer was 
immediately deposited without breaking vacuum. Two 
contacts were made to the film using Eccobond Solder 
59C. Parylene-C was deposited to a thickness of 590 nm, 
conformally encasing the sample on both the top and 
bottom surfaces with a measured thickness uniformity 
better than 10 nm. NiCr films were evaporated onto the 
parylene-C, the bottom surface serving as an absorptive 
broadband antireflection coating and the top surface as 
a gate patterned by shadow masking. The antireflection 
coating was deposited to a sheet resistance of 275 fi/D 
and the gate to 400 fl/D. Two contacts to the gate were 
made using the same silver epoxy. 

Parylene-C optical properties were characterized uti- 
lizing FTIR transmission measurements. The static di- 
electric constant was measured with a capacitance bridge 
and found to be 2.8 for par ylene-C at low temperature 
consistent with reference. |30| The gate provides 2.6 x 10 10 
e/cm 2 per volt. Breakdown voltages as large as 200 V 
were achieved. 

FTIR measurements were performed with a Bomem 
DA8. A far-infrared laser cavity pumped by a CO2 
laser provided the fixed frequency source for the Faraday 
and cyclotron resonance measurements. The polarization 
modulation technique used to measure the Faraday angle 
is detailed elsewhere^! Circular polarized light for the 
cyclotron resonance measurement was generated using 
quartz quarter-waveplates with NiCr antireflection coat- 
ings. A GaAs 2-DEG heterostructure was used to set the 
wavcplatc angle, verify the retardance of the waveplatcs, 
and served to calibrate the Faraday angle measurements 
and the sign of the charge carriers in the cyclotron reso- 
nance measurements. 



edge. The average of the two dispersions is used in the 
analysis and is given by E avg = 1.9fc + 12.6fc 2 + 2, 300/c 6 . 
The TSS Fermi velocity vf = dE/d(Hk), cyclotron mass 
m c = fikp/vp, and carrier density n = k F /(4ir) are used 
in analyses. 



Back surface accumulation layer 

The measured bulk mass near the conduction band 
edge is 0.15 m^S^m The ARPES dispersion gives a 
mass, shown in FigureUJc), of 0.14 too. Both masses are 
expected to increase with Fermi level and remain similar 
in the vicinity of the conduction band edge i 14 ' 15 High in 
the conduction band, the two dispersions become very 
similar so the cyclotron masses are expected to be nearly 
the same. 

The large carrier concentration measured by the zero- 
gate optical measurements in Figure [1] n « 10 13 enr 2 , 
indicates that the Fermi level is high above the conduc- 
tion band edge, dc transport measurements on similar 
films show nearly the same 2-D carrier density over a 
wide range of film thicknesses 

If these carriers were uniformly distributed across the 
60 nm film, the resulting bulk carrier density would 
be rib ~ 2 x 10 18 cm" 3 . This value should be com- 
pared with single crystals that are routinely an order 
of magnitude lower in bulk carrier density^ Assuming 
E F = (hk F ) 2 /2m* and n b = k F /{3ir 2 ) where the band 
edge effective mass is take to be 0.15 mo, the resulting 
Fermi energy above the conduction band edge is 34 meV. 

As mentioned in the main text, dc transport argues 
against this interpretation since similarly grown epitax- 
ial films show the bulk has a very small contribution to 
the conductivity compared to the accumulation layers for 
such thin films ~ The observed large carrier density is 
shown to be associated with a bottom surface accumula- 
tion layer likely due to defects which can result in a very 
high Fermi level. 



II. THOMAS-FERMI SCREENING AND 
OPTICAL MODELS 

Topological surface state dispersion 

Peaks associated with the topological surface state 
in the momentum distribution curves as measured by 
ARPESi&±£ were fit along the T-K and T-M directions 
at many binding energies ranging from the valence band 
to well above the conduction band edge. Fits to the dis- 
persion where energy is measured from the Dirac point 
give E\ T -k = l-7fc + 13.5fc 2 - 4,600fc 6 and E\ T - M = 
2.2fc + 8.2£; 2 + 4, 200fc 6 where E F is in eV and k is in A~\ 
There is little difference between the two dispersions from 
hexagonal warping 3 -^ since our range of Fermi level is 
from the Dirac point to just above the conduction band 



Thomas-Fermi screening model in topological 
insulators 

Knowing the capacitance, a gate transfers a deter- 
mined amount of net charge to or from a topological insu- 
lator. The net charge is distributed between the topolog- 
ical surface state and the bulk. Deriving this proportion 
requires solving the general band bending problem in a 
topological insulator due to screening of potentials^ 2 - 

In the Thomas-Fermi screening model, the potential 
satisfies Poisson's equation which we solve for a semi- 
infinite slab of topological insulator subject to an exter- 
nally applied electric field. Figure [ST] shows band bending 
of the conduction band edge as a function of depth with 
a Fermi level Ep^ defined away from the screened region. 
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Depletion Length 




FIG. SI: Calculated band bending effects of the conduction 
band edge (red, blue, purple) in a semi-infinite slab of topo- 
logical insulator. Ept (dashed green) is the Fermi level deep 
in the bulk, (f> is the potential due to an externally applied 
electric field, A s is the screening length, and x is the distance 
away from the gated surface. Depletion onset occurs when 
ecj> — EFb- For larger cf> (red), a depletion layer forms separat- 
ing the topological surface state from the bulk carriers. Zero 
<j> (blue) defines flat band where the background ion density 
equals the electron density. For negative <j> (purple), the con- 
duction band edge bends downward creating an accumulation 
layer. 



d 2 (j)(x) p(x) 
dx 2 e 

The net electrostatic charge density is related to the 
net ionic background density n and the carrier density 
of electrons n(x) by p{x) = e(no — n(x)) where no — 
^(2mE Fb /h 2 f/ 2 and n(x) = £z(^(E Fb - e<j)(x))) 3 / 2 

The differential equation can be written in the follow- 
ing general form: 



where 5 



= (l-(l-5) 3/2 ) 



f-, and A 2 



A, defines the 



d 2 S 
dy 2 

"ett' y — 

screening length and e 

^ and integrating, the second order differential equation 
is converted into a first order one. For the specific case of 
a semi-infinite slab where Es and 5 are zero deep inside 
the topological insulator, the solution is simply given by: 



100 for Bi 2 Se 3 ^ Defining E s 



Ec 



±(25+ -(1-5) 



5/2 



,1/2 



Where the + (-) sign is for the case of band bending 
down (up), or equivalently 5(0) < (5(0) > 0). Deple- 
tion onset occurs when 5(0) = 1. Either Es or 5 on the 
inside surface at y = is determined given an externally 
applied gate voltage and the uncompensated charge asso- 
ciated with the topological surface state whose dispersion 
is characterized by ARPES. 



Optical Models 

The general optical models which were used in sim- 
ulating the various data sets are described. Approxi- 
mate formulas relating to various optical measurements 
are presented first which emphasize the physical concepts 
with minimal interference from details. 

The normal-incident transmission T for a thin film 
with multiple contributions to the conductivity a rela- 
tive to the transmission of the substrate To is given by 

= | tttv' — 12 where v; = ZnCTi is the surface ad- 

mittance, n su b is the index of refraction of the substrate, 
and Zq is the impedance of free space. 

Fourier transform spectroscopy (FTIR) measurements 
performed in zero magnetic field are described within 
a Drude model defined by <7 Tr = " p3 ^ Z " where u>„„ = 
Z Q ne 2 /m. a xx is the longitudinal conductivity, u) ps is 
the two dimensional plasma frequency, 7 is the scatter- 
ing rate, 10 is the radiation frequency, n is the two dimen- 
sional carrier density, e is the electronic charge, and m is 
the effective mass. 

Cyclotron resonance (CR) transmission measurements 
are performed at fixed frequency with normally inci- 
dent circularly polarized light function of applied 
magnetic field. The surface admittance terms are most 
simply given in the circular polarization basis where 
± _ — u P s/z w } iere u — eB/cm r is the cyclotron 

7 — l(OJ±LJ c ) L/C J 

frequency, m c is the cyclotron mass, and er* is the con- 
ductivity in the circular polarization basis. In the lim- 
iting case yi « n su b + 1, the transmission reduces to 
T ~ „ 2 ^+i\ 2 I 2 %o I • The cyclotron resonances can 
be thought of as a sum of absorptive Lorentzian reso- 
nances centered at B res .i = (cuj / e)m c ^. Differential cy- 
clotron resonance measurements (A-CR) can be thought 
of as a sum of differences of Lorentzians. The polarity 
and value of the B-field resonance determines the sign 
and mass of the carriers. 

The complex Faraday angle is defined as tamOp = 



1 v ~ i 
L t++t- 



where tanflp « 9p, t are the 



transmission Frcsncl coefficients in the circular polariza- 
tion basis, and 2ia xy = <j + — o~ *21 \miQp) is related to 
the induced circular dichroism similar to cyclotron reso- 
nance measurements. The real part, thought of as the an- 
gle of polarization rotation, is related to the reactive part 
of the conductivity. When the scattering rate is much less 
than the cyclotron frequency, lm(0p) (Re(6p)) shows a 
simple Lorentzian resonance (antiresonance) centered at 
B res . Contrary to CR measurements, n and p-type car- 
rier resonances will simultaneously appear in both polar- 
ities of magnetic field but with opposite sign. 

The actual model used to simulate the data is derived 
from a general solution to the transmission Fresnel co- 
efficient for light normally incident on a stack of slabs 
which may have multiple conductingfilms at interfaces. 
The formalism is given in Reference [3, Appendix I. 

Maxwell's equations inside media with appropriate 
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boundary conditions are used to derive two matrices 
which can be viewed as transfer functions for the elec- 
tric field vector .— The 2x2 matrix M [m , n,+i , cr] transfers 
the electric field across an interface consisting of a film(s) 
with a net 2D conductivity a from a slab of index rii to 
a slab of index n 4+ i. The matrix K[rii,di] propagates 
the electric field through a particular slab of thickness di 
with index rij. 

Since the parylene coating and In 2 Se3 capping layer 
are thin (compared to wavelength) non-absorbing dielec- 
tric films in our frequency range, they do not effect the 
optical signals so we leave them out. The Bi2Se3 film 
is idealized as a conducting 2D film on top of an insu- 
lating Bi 2 Se3 dielectric slab that has optical absorptions 
due to phonons. Therefore the first boundary is between 
air and the idealized Bi2Se3 insulating slab, and between 
them is a NiCr film (used as a gate and characterized by 
the dc sheet resistance R g ) and another conductor asso- 
ciated with the Bi2Se3 free carrier response represented 
as (T ± . All contributions from all the free carriers in the 
Bi2Se3 film are lumped into cr*. The second interface be- 
tween the Bi2Se3 insulating slab and sapphire substrate 
is considered to have no free carriers. The last interface 
is sapphire into air with a NiCr film, an anti-reflection 
coating characterized by sheet resistance Rar- 

Using the transfer matrices, the Fresnel transmission 
coefficients are derived for the entire stack of slabs. All 
optical quantities presented in this paper are calculated 
from this expression. 

The cyclotron resonance transmission is given by 
T + [V g ] = t+t+. The FTIR transmission spectra, al- 
though unpolarizcd, is given by the same expression with 
B = and o~ + replaced by the longitudinal polariza- 
tion expression o~ xx . The Faraday expression is given by 

9 F [V g ] = i£=£. 

For the differentially gated measurements, A-CR is 
given by 

AT+ T+[V g + AV/2] - T+[V g - AV/2] 
T^g~ ~ \{T+[V g + Ay/2] + T+[V g - AV/2]) 

Unpolarizcd gated FTIR measurements are reported nor- 
malized to the zero gate transmittance and is given 
by (where B is set to zero) T +ty l i ■ A9p is given by 
A6 F [Vg} = 9 F [Vg + AV/2] - e F [V g - AV/2]. 

There are two prominent phonons in Bi2Se3 over our 
spectral range of interest which can be written in terms 
of the complex index of refraction nsiSe = \[i where 

2 n 2 

e = eo + J2 e i and e i = ^-^-iur- ■ Tlie P nonon param- 

4=1 

eters for epitaxial films are experimentally determined 
to be eo = 25.6, Oi = 615 cm -1 , u\ = 63.1 enr 1 , 
Ti = 1.8 cm" 1 , VL2 = 80 cm" 1 , u>2 = 133 cm" 1 , and 
T2 = 2 cm -1 . Similar values for bulk crystals were pre- 
viously reported.— 

The free carrier response of the Bi2Se3 film is modeled 

3 , 

as the sum of three Drude terms cr* = V — "p°,'/ - — T 



where w Cj i = eB/cm c ,i is the cyclotron frequency, uj ps .i = 
Zoriie 2 1 'm c ,i is the two dimensional plasma frequency, 7, 
is the inverse transport lifetime, m c .i is the cyclotron 
mass, rii is the two dimensional carrier density, and Zo is 
the impedance of free space. 

The parameters for the Drude term labelled BSAL 
(bottom surface accumulation layer) in the main text is 
found to be n(, s = 1 x 10 13 cm" 2 , jb s = 50 cm" 1 , and 
mbs = 0.25 mo where mo is the mass of a free electron. 
For the Drude term labelled TSS (top surface topological 
surface state), rnrss and n-TSS are set equal to the values 
based upon ARPES digitized TSS dispersion such that 
n-TSS = k F /(4ir) and ititss = hkp/vp where k F and 
v F are functions of Ep. Ep is a function of gate volt- 
age determined by the Thomas-Fermi screening model. 
The relation between gate voltage V g and Fermi level Ep 
are severely constrained by the small mass observed in 
the raw data at the highest negative gate voltages. The 
TSS scattering rate jtss is a free parameter which is 
measured. 

The surface admittance associated with the carriers 
labelled IB (interim bulk) is assumed to be much smaller 
than the BSAL contribution in this experiment. In this 
case, the value chosen for the carrier density associated 
with the IB carriers at the Dirac point is insignificant in 
the overall differential optical response. Only the change 
in carrier density is relevant. We therefore assume for 
convenience that nis = at the Dirac point. All changes 
in carrier density of the IB carriers are then completely 
determined by the band bending model. The scattering 
rate 7/5 and cyclotron mass mis are free parameters. 

NiCr films are used as a gate and antircflcction (AR) 
coating since NiCr has a very high scattering rate. There- 
fore, the admittance is frequency independent at THz fre- 
quencies. The optical response is then completely char- 
acterized by the dc sheet resistance where Rar = 275 fi, 
and Rg = 400 f2. Other parameters are the index of re- 
fraction of sapphire n<? a pft — 3-1 and vacuum no = 1, the 
thickness of the sapphire substrate dsaph = -05 cm, and 
the Bi2Sc3 film thickness ds,se = 60 nm. 

The amount of band bending at zero gate voltage as- 
sociated with defects can be conceptualized as caused by 
fixed surface charge. It is a free parameter, but highly 
constrained by the observed small mass at the highest 
negative gate voltage combined with knowing the amount 
of charge transferred from a gate into the surface state, 
and is found to be 2.8 x 10 12 e cm~ 2 for the top surface. 

It is not a priori obvious that three Drude terms should 
suffice in reproducing the data since numerous conductiv- 
ity contributions could be associated with bulk and accu- 
mulated carriers. For the BSAL, there exists many sub- 
band states, and the lowest subbands are likely Rashba 
spliti 16 ' 18 Although the mass of the IB carriers (0.21 ) is 
expected to be between the conduction band edge mass 
(0.15)£22 and the BSAL mass (0.25), it is quite large and 
likely a result of a convolution involving an average over 
bulk carriers and low binding energy BSAL states. How- 
ever, the differential optical model that is insensitive to 



9 



the BSAL parameters impressively reproduces the data 
sets with effectively two Drude terms allowing clean ex- 
traction of the TSS parameters. 



III. GATED OPTICAL DATA ADDENDA AND 
OPTICAL MODEL DISCUSSION 

Gated optical data addenda 

Figures IS2fa.b) report A-CR data that are similar to 
Figure [Ha) except the measurements were performed at 
different frequencies. Figurcs lS2T d.e) are the modeled op- 
tical responses that use the same parameters as reported 
in Figure [H 

Differential Fourier transform spectroscopy data (A- 
FTIR) and modeled response are reported in Figures 
IS2T c) and (f). The oscillations present in the data are 
from Fabry-Perot interference occurring in the substrate 
due to an imperfect antireflcction coating. The increase 
in transmission with negative gate voltage is evidence of 
dominant n-type carriers. If a single Drude n-type carrier 
existed and a gate changed only the carrier density, then 
the data set would be fully symmetric about zero gate 
bias. The asymmetry is an indication of surface carrier 
contributions changing scattering rate and/or mass. 

A 7 mm spot size was used for the FTIR measurements. 
Such a large spot size relative to the sample allows the 
maximum of throughput power required to measure to 
low frequencies. However, some leakage light occurred 
through ungated parts of the film due to the shadow 
masking procedure. The measured A-FTIR signals are 
therefore reduced. Although the photometries differ from 
the model, the qualitative behavior is reproduced. The 
apparent asymmetric response observed in the transmis- 
sion associated with the phonon resonance in FigurcsHJc) 
and lS2f c) is reproduced by a simple Lorentzian dielectric 
response function and requires no other exotic explana- 
tion. 



Intuitive description of Figure 2(d), A-CR model 

The modeled optical signals in relation to the Drude 
parameters can be qualitatively understood most easily 
by considering the A-CR model curves of Figure 2(d). A 
CR is a Lorentzian-like absorption where n defines the 
weight, 7 defines the width, and m c defines the center of 
the resonance via B res = (cuj j e)m c . When a gate voltage 
is applied, these parameters can change. A single A-CR 
curve is the difference between pairs of Lorentzians, a 
pair for the TSS and IB contributions. It is convenient 
to ignore the BSAL completely since it does not change 
with gate voltage and therefore does not contribute to 
the differential response. Also note that the surface ad- 
mittance associated with the BSAL contribution is much 
larger than the IB and TSS contribution, so the overall 



modeled differential optical signals arc insensitive to the 
precise values of the BSAL parameters. 

Near zero gate voltage, the mass and scattering rate 
of the TSS and IB carriers do not change much with 
gate voltage. Only the weights of the resonances change. 
Since jtss >> Iib and the change in carrier density 
Anrss ~ An/s as shown in Figure [3j the resulting line 
shape appears Lorentzian centered on the mass of the 
IB contribution. As the gate is increased and the Fermi 
level moves higher, the dominant response is the IB con- 
tribution that has a growing scattering rate and carrier 
density. Taking the difference between two Lorentzians 
that have slightly different widths (7) will cause a sup- 
pression of the peak at the resonant B-field. A slight si- 
multaneous increase in mass shifts one of the resonances 
in B-ficld relative to the other causing a slight skew in the 
difference curve. In this case, the part of the differential 
cyclotron resonance occurring at low B-fields arc slightly 
more suppressed than the higher fields due to the chang- 
ing IB mass. The gradual increase of the TSS scattering 
rate which is relatively large causes a much wider gen- 
eral suppression, most notably giving rise to the subtle 
suppression of the tails at large negative B-fields. 

As one gates from zero voltage toward the Dirac point, 
only the carrier density of the IB contribution changes. 
Therefore, the A-CR of the IB carriers maintain an ap- 
proximately constant Lorentzian— like line shape. The 
TSS response is superimposed on the IB background. A 
rapid decrease of the mass in the TSS contribution in 
combination with the rapidly diminishing spectral weight 
give rise to the peak-dip-hump structure. The change in 
scattering rate amplify these features. 

At the highest negative gate voltages, the A-CR peak 
for the TSS becomes just the CR peak of the higher Fermi 
level since very little relative spectral weight exists in the 
lower Fermi level peak. This is most vividly illustrated 
by considering the hypothetical scenario where the lowest 
Fermi level is actually at the Dirac point where there are 
no carriers. The difference curve then would be exactly 
the TSS CR curve at the higher Fermi level, superim- 
posed on the IB contribution. Therefore, the peak-dip- 
hump structure evolves to a single peak structure at the 
two highest negative gate voltages. 

An understanding of all the A-CR curves can be sum- 
marized. The IB contribution is the dominant response 
at positive gate voltages. At negative voltages, the IB 
contribution to A-CR is roughly constant since only n is 
changing. The peak-dip-hump structure that develops at 
large negative gate voltages is due to the rapidly chang- 
ing mass amplified by the changing scattering rate of the 
TSS. At the highest negative gate voltages, a single TSS 
cyclotron resonance is superimposed on a constant A-CR 
IB Lorentzian-like peak. 
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FIG. S2: Gated optical transmission measurements addenda (a-b) 10K A-CR measurements are the same as described 
in the caption of Figure [2] except measured at different frequencies labeled on each graph (c) 6 K zero magnetic field FTIR 
transmission measured at gate voltage V g normalized to the zero gate value, (d-f) the corresponding modeled optical responses 
using the same parameters reported in Figure [3] that reproduce the data of Figures Oa-c) and Figures [TJc,d). 



IV. TOPOLOGICAL INTERFACE STATES 

The interface between two insulating materials with 
band gaps inverted with respect to each other can sup- 
port topologically protected interface states^^ For the 
ARPES measurements on bulk topological insulating ma- 
terials the vacuum takes the part of the topologically 
trivial insulator. The more general case has been solved 
within the two band model using k ■ p perturbation the- 
ory in reference HH This theory shows that the Dirac 
spectrum is modified as a function of the gap of the 
trivial insulator and any potential step at the interface. 
Thus for a trivial insulator with a gap of 1.2 eV in con- 



tact with a topological insulator with a bulk gap of 300 
meV in the presence of a potential step of 0.3 eV the 
Dirac point of the topological surface state shifts by 60 
meV. This potential step can result from differences in 
the work function of the two materials or the presence 
of a dipole layer at the interface. While not enough is 
known about the In2Sc3/Bi2Se3 interface to predict the 
potential step^~— this calculation demonstrates that a 
shift of the right magnitude is expected to occur with 
reasonable values of a potential step. A more accurate 
theory together with experimental determination of the 
potential step is needed to confirm this result. 



